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Editorial 


As Zoological Research enters its 33rd year of publication, let us reflect back on the last three decades of dramatic 
change to China and the scientific community. In the latter half of the 20th century, scientific research was almost 
exclusively published in the United States or Europe. We started Zoological Research in 1980 with the goal of creating 
a strong Chinese-language journal to support regional and national research and push towards developing an innovative 
and collaborative Chinese scientific community. Over the next two decades, we worked diligently to make Zoological 
Research in to one of China’s top life science journals, soliciting manuscripts from across the country. By the mid 1990s, 
however, we realized these efforts were not enough. Chinese researchers and research on China had begun to play a 
much larger role in the international scientific community. Accordingly, we decided to open up submissions to include 
articles published in English by both Chinese and foreign scientists. In 1996, we published our first set of English 
articles, all by Chinese authors, detailing regional findings: a new genus of dragonfly in Tibet, the nesting patterns of 
ants, and a study of snake karyotypes from Yunnan province. 

After opening submissions to both English and Chinese articles, Zoological Research began to gain international 
traction, publishing articles in both languages on a variety of subjects, areas, and from researchers across the globe. 
English-language submissions to the Zoological Research steadily increased, mirroring China’s rise to become a major 
player in the global scientific community as the no. 2 global producer of research, accounting for 14% of all research in 
2010, publishing over 200,000 articles annually. Our journal has seen similar trends, embracing a strong spirit of 
international outreach and cooperation by annually publishing a third of our 100+ articles in English. We have likewise 
published two special issues on primates and animal models of human diseases and are now drawing submissions 
world-wide in preparation for a third. In 2010, we reached a critical milestone when Zoological Research became 
indexed by PubMed/Medline, in addition to Biological Abstracts, Zoological Record, and Chemical Abstracts, among 
others. Likewise, today nearly a third of our 66 member editorial review board is based outside China. 

In the growing spirit of international cooperation that has so deeply enriched the scholarship of our journal, we are 
pleased to announce that starting this month, Zoological Research is a full dual-language publication, concurrently 
publishing Chinese and English issues. We moved to this new format in order to take advantages of our journals two 
great strengths: a strong regional focus of the finest Chinese scholarship and a growing, internationally focused array of 
English-language research submitted by authors around the globe. We hope this new format maintains our distinction as 
one of China’s leading zoological journals while also embracing a new international outlook. To aid in this endeavor, 
our sponsoring institution, the Kunming Institute of Zoology, has hired an in-house language editor from the United 
States who will be offering complimentary English language editing on incoming research and working closely with our 
submitting authors on refining their manuscripts for publication in the new Zoological Research. 

At this time, we would also like to make a special invitation to our readers who are interested in serving as single- 
issue Guest Editors for upcoming issues. Researchers interested in organizing special themed issues to explore a 
particular topic are heartily invited to contact our editor-in-chief, Zhang Ya-ping for more details. As always, thank you 
for your constant support of our journal and our work. Please continue sending your prospective manuscripts and 
following our journal as we continue in our efforts to make Zoological Research a globally known and respected forum 
for zoological research. 


Best, 


Yaping ZHANG 
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Genetic diversity and differentiation of masu salmon (Oncorhynchus 


masou masou) between and within cultured populations 


inferred from microsatellite DNA analysis 


Zhiying JIA’, Yuyong ZHANG’, Shugiang CHEN" ?, Lianyu SHI'” 


1. Heilongjiang River Fisheries Research Institute, Chinese Academy of Fishery Sciences, Harbin 150070, China; 
2. College of Fisheries and Life Science, Shanghai Ocean University, Shanghai 201306, China 


Abstract: Masu salmon, Oncorhynchus masou masou, is one of the most valuable fishery species that has been introduced to China, 
though to date no studies on the genetic diversity and genetic relationship among hatchery populations has been performed with 
molecular markers. We undertook such a study and sampled 120 individuals from three hatchery stocks and analyzed 20 
microsatellite loci. All loci were polymorphic and a total of 91 alleles were detected. A relatively low level of genetic diversity was 
revealed with effective number of allele of 3.1094, 3.3299 and 3.1894 and expected heterozygosity of 0.6600, 0.6648 and 0.6638 in 
the three stocks, respectively. Deviations from Hardy-Weinberg equilibrium were found due to heterozygote deficit. Accordingly, 
evidence of genetic bottlenecks were found in the three stocks. An individual assignment test demonstrated that 8596 of individuals 
were correctly assigned into their original stocks. Pairwise Fst revealed that significant differentiation occurred between these three 
stocks. The results of the study indicated that disequilibrium of genetic structure and differentiation has occurred in all three stocks. 
This information collectively provides a basis for measures to avoid of loss of genetic diversity and introgression in Chinese 


aquaculture. 


Keywords: Oncorhynchus masou masou; Microsatellite; Genetic diversity; Genetic differentiation 


Masu salmon, O. masou masou, were widely 
distributed in Japanese mountain streams (Kano et al, 
2010) before being introduced to China in the 1990s 
(Wang, 1998). The fish is now one of the most valuable 
fishery species in cold-water aquaculture. At present, 
commercial catches of O. masou masou in Japan have 
decreased and hatchery programs have been established 
to release fry to its habitat (Yu et al, 2010). In China, 
because of its good taste and flesh quality, the fish has 
been an important aquaculture species (Wang, 1998). 
However, successive artificial selection and serious 
stress from inbreeding have led to the remarkable 
degradation in the fish: early sex maturity, smaller size 
and silvered body color (Zhang et al, 2011). A second 
batch of O. masou masou was introduced into China in 
2008 to counter these changes. 

High genetic diversity is important for the population 
conservation because it is closely related to population’s 
fitness and adaptability to environment as well as 
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selection (Filgueira et al, 2010). Loss of genetic variation 
among hatchery stocks is a common phenomenon caused 
by genetic drift in many species, including mrigal (Aung 
et al, 2010), Japanese flounder (Sekino et al, 2002) and 
common carp (Kohlmann et al, 2005). At present, there 
are three O. masou masou hatchery stocks in China, but 
their genetic diversity and relationship have not been 
reported despite clear indications of genetic 
abnormalities that prompted the second importation of 
stock in 2008. To avoid a repeat of the earlier incidence, 
information on genetic variation of these hatchery stocks 
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is urgently needed to develop a coherent management 
and conservation plan for this species. 

Microsatellites DNA markers have been proven to 
be useful markers in researching genetic diversity and 
population structure in fish such as Japanese flounder 
(Sekino et al, 2002) and common carp (Kohlmann et al, 
2005). We aim to investigate the genetic diversity and 
population relationships of the Chinese three stocks of O. 
masou masou using 20 pairs of microsatellite primers 
developed for rainbow trout to assess the genetic 
diversity and their relationship within and among stocks 
of O. masou masou. 


MATERIALS AND METHODS 


Sample collection and history of stocks studied 

Two batches of fertilized eggs of O. masou masou 
were introduced into China, the first in 1996 and 
the second 2008. Both were kept in  Bohai 
Experimental Station. In 1998, parts of O. masou 
masou individuals were taken to Beijing Shuntong 
Rainbow Trout Culture Center and regarded as a separate 
stock population after several generations of artificial 
propagation, resulting in three stocks. We named the 
three stocks we used BO1996, Beijing, and BO2008. In 
2010, we collected 40 samples from each of the 
presumptive stocks. And the samples of the BO2008 
stock were collected directly from the introduced 
population. Fin-clip tissues were preserved in 75% 
ethanol (Jia et al, 2008) for BO1996 and Beijing stocks. 
ВО2008 stock samples were collected as fry just before 
beginning to float. 


DNA extraction and microsatellite DNA analysis 

Genome DNA was extracted using the standard 
proteinase К, phenol-chloroform-isoamyl procedure 
(Wei et al, 2001). 20 microsatellite loci (Tablel) 
developed for other salmon were used in the study. The 
OMM3005 primer was from NCBI (http://www. 
ncbi.nlm.nih.gov/nuccore/20142355) and others were 
from the articles of Rexroad et al (2002a,b;2005). Primer 
conditions, mainly annealing temperatures, were 
optimized for amplification. PCR was performed using a 
reaction volume of 15 uL containing 50 ng of genomic 
DNA, 0.2 umol/L of each dNTP (TaKaRa, China), 
0.1 mmol/L of each primer, 1.5 mmol/L MgCl, 1х 
PCR buffer (TaKaRa, China) and 1U Taq DNA 
polymerase. The thermal cycling (Biometra) began with 
3 min at 94 °C, followed by 27 cycles of denaturing at 94 
°C for 30 s, annealing temperature from 50 °C to 62 °C 
depending on different primers for 30 s, elongation for 
30 $ at 72 °C, followed by a final extension for 5 min at 
72 °С. РСЕ products were separated іп 10% 
polyacrylamide gels. 
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Statistical analyses 

The effective number of alleles per locus (Ne), 
number of detected alleles (Na), observed heterozygosity 
(Ho), expected heterozygosity (He) and allele frequency 
were calculated using POPGENE 1.32(Yeh et al, 2000). 
To test the difference between populations of Ne and He, 
non-parametric t-tests (Sokal & Rohlf, 1995) were 
performed using Microsoft Excel. 

Deviations from Hardy-Weinberg equilibrium 
(HWE) were estimated using GENEPOP 4.0 (Raymond 
& Rousset, 1995). Exact P-values were evaluated with 
Markov chain method (dememorization=10000; 
batches=100; iterations per batch=5 000). Additionally, 
deviations from HWE in the direction of heterozygote 
excess or deficit were also carried out. 

Analysis of molecular variance (AMOVA), 
implemented in GenAIEx (Peakall & Smouse, 2006), 
was used to partition total genetic variation hierarchically 
among stocks. 

Genetic distance between populations was estimated 
by the PHYLIP software package, version 3.69 
(Felsenstein, 1993). Bootstrap re-sampling (1 000 
replicates) was performed to test the reliability of 
dendrogram implemented in PHYLIP software. 

STRUCTURE 2.2 (Pritchard et al, 2000) was used 
to determine the most likely number of genetic clusters 
(К=1-5). An assignment test was also performed in 
order to determine the extent to which individuals could 
be correctly assigned to their population of origin. In the 
analysis, the simulation was used assuming admixture 
and correlated allele frequencies between populations 
with a 10000 replications burn-in period and 100 000 
MCMC replicates. 

Evidence of population bottlenecks were tested with 
Bottleneck 1.2.02 (Cornuet & Luikart, 1997) and the 
Wilcoxon sign-rank test was used to determine if a 
population exhibits a significant number of loci with 
heterozygote excess with a two-phased model of 
mutation (TPM, as recommended in the program). The 
mode-shift indicator was also employed to discriminate 
bottlenecked populations from stable populations. 


RESULTS 


PCR results of 20 pairs of microsattelite primers 
According to the clarity and size of banding 
patterns, 20 pairs of primers showed good images and all 
were polymorphic in the sampled populations (Figure 1). 


Genetic diversity and linkage disequilibrium 
Summary statistics showing allelic variation within 
loci and population are presented in Table 1. A total of 91 
alleles were detected in the 20 loci among 120 
individuals with an average of 4.55, ranging from three 
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Figure 1 Polyacrylamide gel results at the ОММ5017 locus for three stocks of masu salmon 


a, b, c, d and e means different alleles. 


alleles at loci AF352740, AF352758 and ОММ5020 to 7 
alleles at locus AF352763. All loci were polymorphic in 
all samples from three stocks. Ho and He within stocks 
ranged from 0.2973 (BO2008) to 0.9714 (BO1996) 
and from 0.3225 (Beijing) to 0.8290 (Beijing), 
respectively. 

Exact tests for genotypic linkage disequilibrium 
showed significant deviations in the stocks of BO1996 
(OMM1036 and OMM1035, OMM1036 and AF352744, 
AF352738 and AF352750 pairs) and Beijing (OMM1045 
and AF352763, AF352763 and af352746, AF352750 and 
AF352744 pairs) at the 5% level. But disequilibria were 
not consistently observed in the same locus pairs in each 
population. 


Genetic structure with and among populations 

The effective number of allele per population was 
3.1094, 3.3299 and 3.1894 in BO1996, BO2008 and 
Beijing respectively. The mean expected heterozygosity 
per population was 0.6172, 0.6171 and 0.6299 in 
BO1996, BO2008 and Beijing respectively. No 
significant difference of genetic diversity was found 
between stocks. 

Deviations from HWE were found in all three 
populations (Table 1). Heterozygote deficit was observed 
in BO2008 stock while for the other two, however, 
heterozygote deficits existed at only few loci and did not 
achieve population level. The ВО1996 population 
exhibited the highest number of loci (11) deviated from 
HWE and then BO2008 (10) and Beijing (9) stock 
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populations. Two loci (ОММ5017 and ОММ5010) 
deviated from HWE in all populations but a heterozygote 
deficiency existed in different population (ОММ5017 for 
BO2008 and Beijing, ОММ5010 Юг BO2008 at 5% 
level). A significant bottleneck was found for three 
stocks (P [heterozygote excess] = 0.00000, 0.00001 and 
0.00136 for BO1996, BO2008 and Beijing, respectively 
and the shifted model). 

For the real number of clusters, once the cluster (K) 
is reached, L(K) at larger Ks plateaus or continues 
increasing slightly. In our analysis, L(K) at larger Ks 
plateaus when K=3, which is consistent with the groups 
as individual sampled (Figure 2). And the 4K value that 
could give a clear peak at the true value of K according 
to procedures of Evanno et al (2005) was in Figure 3. 
The self-assignment test resulted in correct assignment of 
102 of the 120 individuals in the three sampled groups, 
corresponding to 85.0% (Figure 4). The number of 
misclassified individuals in the three groups was four, 
seven, and seven with the ratios of 3.33%, 5.83%, and, 
5.83% respectively. 

The locality of origin for each individual is 
indicated in the X-axis. 


Genetic relationships 

The global test for population differentiation was 
significant with Fst value of 0.055 (P=0.01). Pairwise 
Fst estimates and level of significance for all population 
pairs are presented in Table 3. All Fst values were 
significant (all P=0.01) and the highest Fst was found 
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Table 1 Genetic variability of three O. masou masou stocks at 20 microsatellite loci 





Locus Repeat Population P Ne Ho He 
РР ВО1996 0.0013* 3.9801 0.7500 0.7582 
(GT)17 BO2008 0.0077* 3.8186 0.6250 0.7475 
BEIJING 0.0317* 2.4078 0.6000 0.5921 
BO1996 0.0000* 3.9851 0.9000 0.7585 
OMM5010 (GT)25 BO2008 0.0000* 4.3185 0.5250 0.7782 
BEIJING 0.0017* 3.8415 0.8250 0.7491 
BO1996 0.0988 3.3058 0.8250 0.7063 
Я (GA)44 BO2008 0.0146* 2.4298 0.7500 0.5959 
BEIJING 0.4824 2.7997 0.7250 0.6509 
BO1996 0.0084* 2.9331 0.6750 0.6674 
OMM5008 (GT)19 BO2008 0.2133 2.4600 0.4865 0.6016 
BEIJING 0.0058* 2.8855 0.8500 0.6617 
BO1996 0.2681 3.0149 0.8205 0.6770 
DMBDUD (CA)18 BO2008 0.0003* 2.5499 0.8462 0.6157 
BEIJING 0.0475* 3.3934 0.8250 0.7142 
BO1996 0.0174* 3.7559 0.4750 0.7430 
OMM1045 (TCTA)12 BO2008 0.0003* 2.3145 0.2973 0.5757 
BEIJING 0.2520 3.0274 0.5500 0.6782 
BO1996 0.0111* 3.7736 0.7000 0.7443 
S (TCCA)14 BO2008 0.0000* 5.2605 0.5526 0.8207 
BEIJING 0.1396 4.5403 0.6154 0.7899 
BO1996 0.9844 3.1232 0.7692 0.6886 
MP8 (TATC)18 BO2008 0.5301 4.8784 0.7368 0.8056 
BEIJING 0.1108 3.5242 0.8750 0.7253 
BO1996 0.0620 4.5134 0.6923 0.7885 
eee (TAGA)31 BO2008 0.0508 5.1200 0.6500 0.8149 
BEIJING 0.0000* 4.3478 0.9500 0.7797 
BO1996 0.5363 3.5455 0.6923 0.7273 
OMM1012 (CA)14 BO2008 0.0000* 3.9661 0.4872 0.7576 
BEIJING 0.0388* 3.3465 0.5789 0.7105 
BO1996 0.3406 1.8827 0.6053 0.4751 
OMM1009 (CT)6N2(CT)4 BO2008 0.1241 2.2467 0.7179 0.5621 
BEIJING 0.3765 2.2651 0.6667 0.5658 
BO1996 0.0424* 3.5794 0.6250 0.7297 
OMM3005 (CA)37 BO2008 0.3607 4.3770 0.6667 0.7816 
BEIJING 0.0479* 4.5007 0.8500 0.7877 
BO1996 0.7945 1.7847 0.4250 0.4453 
A332738 (ТАСА)15 BO2008 0.3727 2.5457 0.7500 0.6149 
BEIJING 0.3994 1.9606 0.5526 0.4965 
BO1996 0.1491 4.3986 0.7714 0.7839 
eee (GATA)26 BO2008 0.0579 6.3908 0.7436 0.8545 
BEIJING 0.0802 5.4688 0.9714 0.8290 
BO1996 0.1019 2.3529 0.5250 0.5823 
AF392750 (GATA)14 BO2008 0.0926 2.4317 0.6923 0.5964 
BEIJING 0.0726 2.4634 0.7250 0.6016 
BO1996 0.0115* 2.8594 0.5000 0.6589 
AF392146 (TATC)7 BO2008 0.0070* 2.9304 0.4750 0.6671 
BEIJING 0.1266 3.3508 0.6000 0.7104 
BO1996 0.0112* 2.8889 0.4615 0.6623 
AP332144 (GATA)19 BO2008 0.0000* 2.7112 0.3077 0.6394 
BEIJING 0.0558 3.0239 0.5128 0.6780 
BO1996 0.0000* 2.6644 0.8500 0.6326 
aerials (AC)14 BO2008 0.2775 2.2028 0.5897 0.5531 
BEIJING 0.0000* 2.8803 0.7750 0.6611 
BO1996 0.0114* 1.9987 0.7179 0.5062 
Bese (GATA)13 BO2008 1.0000 1.8376 0.4865 0.4621 
BEIJING 0.4266 1.4672 0.3000 0.3225 
BO1996 0.9234 1.8470 0.5128 0.4645 
OMM5925 (CA)17 BO2008 0.0457* 1.8075 0.6154 0.4525 
BEIJING 0.0000* 2.2923 0.9500 0.5709 
BO1996 3.1094 0.6647 0.6600 
Mean я BO2008 3.3299 0.6001 0.6648 
BEIJING о -=== 3.1894 0.7149 0.6638 


P: estimate of probability of deviation from Hardy-Weinberg equilibrium; *: significant departure from Hardy-Weinberg equilibrium (P<0.05); Ne: the effective 
number of alleles; Ho: observed heterozygosity; He: expected heterozygosity. 
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between the BO1996 and BO2008 populations. Global 
subdivision with the use of AMOVA revealed that the 
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percentage of variation among populations was 5%, 4% 
among individuals, and 91% within individuals. This 
indicates that the total variation of microsatellite loci was 
due to variation among individuals, rather than among 
stocks. 
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Figure 3 AK value estimated at the К ranged from 1 to 5 





Beijing 


Figure 4 Proportional membership (Q) of each individual іп the three clusters 


Table 3 Pairwise Population Fst Values at 20 
microsatellite loci 


Pop ID BO1996 BO2008 Beijing 
BO1996 ажа 0.01 0.01 
BO2008 0.062 aee 0.01 
Beijing 0.050 0.052 ажа 


Fst Values below diagonal. Probability values based on 99 permutations are 
shown above diagonal. 


Pairwise genetic distances among hatchery 
populations (Table 4) ranged from 0.1282 (В01996- 
Beijing) to 0.1349 (BO1996-BO2008), which was 
consistent with Fst results. 


Table 4 Nei's genetic distance (below diagonal) 


Pop ID BO1996 BO2008 Beijing 

BO1996 = 

ВО2008 0.1671 жәж 000002. 

Beijing 0.1282 0.1349 жж 
DISCUSSIONS 


Our study provides useful insights into the genetic 
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diversity of O. masou masou and the data indicates 
relative low levels of genetic diversity but significant 
population differentiation the three O. masou masou 
stocks. 


Within-population diversity and population structure 

Low levels of allele diversity among O. masou 
masou were suggested by the 20 microsatellite loci. The 
number of alleles ranged from 3.1094 to 3.3299 per locus, 
compared to the average value from 12 species of 
freshwater fish (4—9.1) (De Woody & Avise, 2000) and 
from O. masou masou populations in Japanese and 
Korean populations (4=6.42—9.71) (Yu et al, 2010). 
However, the observed and expected heterozygosity 
(Но=0.6001—0.7149 , Не-0.6600-0.6648% was 
moderate, slightly higher than the average of freshwater 
fish (H=0.54 averaged of 13 species) (De Woody & 
Avise, 2000) and lower than О. тазои masou 
(Но=0.698—0.817, Не-0.722-0.814) (Yu et al, 2010). 
For three stocks, samples were introduced as fry and they 
may have comprised a few full-sib families rather than 
representatives of their original populations (Na-Nakorn 
et al, 2004), which resulted in the low level of alleles. 
We also used primers developed for other salmon. Cross- 
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amplified primers and the loss of low-frequency alleles 
after several generations of artificial propagation may 
result in low genetic variation. BO1996 and Beijing 
population showed even lower genetic diversity 
compared to BO2008 in terms of number of alleles and 
heterozygosity. ВО1996 and Beijing populations 
experienced a longer period of being cultured since 
introduction into China as compared to the BO2008 
population. This time difference might lead to a loss of 
diversity as a result of inbreeding and genetic drift. 

In general, a small number of alleles are a signature 
of population bottleneck (Norris et al, 1999) due to allele 
loss. The present results indicate that significant genetic 
bottlenecks were detected with the shifted model in all 
populations, which implies that a dramatic loss of alleles 
happened in our samples. In the next breeding or 
hatcheries, а sufficient number of individuals are 
necessary to avoid serious bottleneck. 

Deviations from HWE were found іп three 
populations due to heterozygote deficiencies, likely caused 
by null alleles, a mixture of independent populations, 
non-random mating, or artificial or natural selection. In 
this study, half of loci deviated from HWE in the 
sampled populations and heterozygote deficits did not 
occur in all the loci and populations, which did not 
reflect features of null allele effect but likely of mixture 
of independent populations (for BO2008 stock), non- 
random mating, or artificial or natural selection (for 
BO1996 and Beijing stock). 


Genetic differentiation between stock populations 

Genetic differentiation was observed between/ 
among different stocks. In all inter-stock tests, deep 
differentiation was found between BO1996 and BO2008 
and then Beijing and BO2008, which was consistent with 
the different stocks’ origin. 

Significant population differentiation was observed 
between BO2008 and the other two. As a newly 
introduced stock, BO2008 might be genetically different 
from original populations. It is interesting that two other 
stocks, BO1996 and Beijing, are also differentiated 
significantly from each other, as Beijing samples were 
recruited from BO1996 stocks. This could be due to 
genetic drift phenomenon such as non-random mating, 
inbreeding or selection that occurred in hatcheries. 
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Identifying individuals to their origin is important 
for breeding programs to improve desirable quantitative 
traits and to conserve gene resources to avoid vast gene 
introgression during broodstock management. In our 
study, the assignment accuracy of individuals to their 
stock of origin was more than 80%. Cornuet et al (1999) 
calculated the accuracy of assignment of an individual to 
the population of its origin according to the number of 
loci used and the number of individuals sampled and the 
degree of population differentiation and pointed out that 
the best combination was 20-30 loci and 8-12 individuals 
per population when the value of Fst was more than 0.05. 
In our analysis with Fst20.05, 20 loci and 40 individuals 
were used, which is sufficient for the assignment. The 
high percentage of correct assignment (8596) in our 
assignment test also demonstrated that the individual 
assignment tests performed well with the level of 
differentiation. 


Implications for the management 

Although the results of moderate heterozygosity 
were revealed, the presence of bottleneck and relatively 
low levels of allele diversity implied that domestication 
had greatly reduced the genetic diversity of O. masou 
masou in China. Sufficient genetic diversity is crucial in 
ensuring the long-term adaptability for domesticated 
stocks of O. masou masou. Consequently, a large number 
of individuals to avoid inbreeding will be necessary for 
the management and conservation of this species. 
Furthermore, we recommend genetic monitoring be 
performed regularly on these and any new stocks in order 
to avoid a loss of genetic diversity. 

Genetic differentiation revealed by Fst was 
observed between/among different stocks, which was 
supported by the real K estimated from STRUCTURE 
analysis. Therefore, genetic dissimilarity existed in the 
different stock of O. masou masou in China and as such, 
it is necessary to maintain the genetic integrity of each 
stock to avoid gene introgression between the genetic 
groups as revealed in the dendrogram. 
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Abstract: The FTZ-F1 genes encode orphan receptors of the nuclear receptor superfamily and in mammals have been found to play 
important roles in the proper development of the adrenal-gonadal axis and sex-determination. We isolated the homologue of FTZ-F1 
in genetically improved farmed tilapia (gfFTZ-F1). The full-length cDNA was isolated from the ovary, which included an open 
reading frame encoding a predicted protein of 486 amino acids. Sequence, tissue distribution and phylogenic analysis of the FTZ-F1 
showed that the gfFTZ-F1 belonged to SF-1/Ad4BP group and that gfFTZ-F1 transcripts were only expressed in the gonads and 
kidney but not in other tissues. Likewise our data suggests that the gfFTZ-F1 gene may share similar functions with other fish and 


mammalian counterparts, though further study is needed to make any definitive conclusions. 


Keywords: GIFT tilapia; FTZ-F1; Cloning; Expression 


Fushi tarazu factor-1(FTZ-F1) is a member of the 
nuclear receptor superfamily originally described as a 
regulator of the Drosophila homeobox segmentation 
gene fushi tarazu (Lavorgna et al, 1991). FTZ-Fl 
homologues have since been recognized in numerous 
species and broadly divided into two subgroups of 
related genes with separate functions and expression 
patterns among higher vertebrates: LRH/FTF (liver 
receptor hormone/a-fetoprotein transcription factor) and 
SF-1/Ad4BP (steroidogenic factor-l/adrenal 4 binding 
protein) (Ellinger-Ziegelbauer et al, 1994). The 
LRH/FTF genes are expressed mainly in the liver and are 
related to the regulation of the a-fetoprotein gene 
(Galameau et al, 1996) and the cholesterol metabolism 
(Nitta et al, 1999) while the SF-1/Ad4BP genes are 
involved in reproductive functions, regulating the 
transcription of many P450 enzymes acting in the 
steroidogenic pathways as well as by modulating the 
development of the hypothalamic pituitary adrenal and 
gonadal axis and the sex differentiation (Hammer & 
Ihgraham, 1999; Parker & Schimmer, 1997). The SF- 
1/Ad4BP genes also regulate the Mullerian inhibiting 
substance (MIS) transcription and promotion of testis 
development (Giuili et al, 1997; Shen et al, 1994). 

In fish, several FTZ-F1 genes have been identified 
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but some forms of FTZ-F1 homologues were not able to 
be classified in to either group (Von Hofsten et al, 2001), 
though the teleost FTZ-F1 may be a regulator in tissue 
differentiation and a factor in sexual maturation. We 
decided to study the gene in genetically improved farmed 
tilapia (GIFT), Oreochromis niloticus, a commercially 
important cultured freshwater fish in which males grow 
faster than females by 40%—50% (Holden & Reed, 1972). 
Regarding the potential role of FTZ-1 homologues in the 
gonadal development and sex differentiation, it would be 
of significance to study the FTZ-F1 homologue in the 
GIFT tilapia. In this study, a homologue cDNA of FTZ- 
F1 was isolated from the ovary of adult GIFT tilapia in 
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order to examine the expression pattern. 


MATERIALS AND METHODS 


Materials 

Adult GIFT tilapia were maintained іп a 
recirculating fresh water tanks at 26+1 'C at the 
Experimental Base of Guangxi Fishery Institute, 
Nanning City, Guangxi Province, China. 

For the cloning and measurement of FTZ-F1 mRNA 
in the GIFT tilapia, the gonads, brain, heart, liver, 
kidneys, spleen, intestinal, gills, muscle from adult 
female and male fish were quickly collected and snap 
frozen immediately in liquid nitrogen and stored at 
—80 °С until use. 





RNA extraction and cDNA cloning 

Total RNA extraction and reverse transcription were 
carried out as described previously (Cao et al, 2010). A 
pair of degenerate primers (Pl: 5'-CGGAACAAGTTC 
GGCccnatnta-3' and P2: 5'-CCACGATGGAGAACAAC 
Gtytgrt-3') were designed according the conserved 
sequences of the FTZ-F1 gene in other teleosts and used 
for isolation of а FTZ-Fl cDNA fragment of 
approximately 600 bp from the GIFT tilapia ovary. 
Amplification was performed at 94 °С for 2 min, 30 
cycles of 94 °С for 30 s, 58 °С for 30 s and 72 ‘С for 1 
min, and 72 °С for 10 min using a BIO-Rad gradient 
PCR thermal cycler. Then, 5'- and 3'-RACE were 
performed to obtain the 5' and 3' cDNA ends of FTZ-F1 
using the 3'-Full RACE Core Set and 5'-Full RACE Kit 
(TaKaRa) according to the manufacturer’s instructions. 
Two gene-specific primers were used for RACE: P3: 5'- 
GCTCCCCTGAGACTCTTACA-3'; P4: 5'- 
ATGCCGAAGTGAC CGTGGAA -3'. 

The amplified products were separated and purified 
with agarose gel DNA purification Kit (TaKaRa). The 
purified fragments were then ligated into pMD 18-T 
Vector (TaKaRa), propagated in Е. coli DH5a. The 
recombinant plasmids were sequenced by Peking 
Audiocodes Biotech Co., Ltd.. 


Sequence analysis and alignment 

The alignment of the FTZ-F1 protein's amino acid 
sequence was performed using the Clustal W algorithm 
1.7 (European Bioinformatic Institute, Hinxton, United 
Kingdom) using default settings (Gap opening penalty 10, 
gap extension penalty 0.05, gap distance 8). A 
phylogenic tree was generated with the neighbor-joining 
method using Mega 4.0. The sequences used Юг 
comparison and phylogenic analysis and their GenBank 
accession numbers were as follows: Acanthopagrus 
schlegelii (sbFTZ-Fla: AAS75791), Acanthopagrus 
schlegelii (sbFTZ-Flb: AAS75792), Bos Taurus 
(btAd4BP: BAA02764), Cynoglossus semilaevis (tsFTZ-F1: 
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ABQ41307), Clarias gariepinus (acFTZ-F la: AAG49004), 
Clarias gariepinus (асЕТ7-Е1Ъ: AAG49005), Carassius 
auratus (Gftz-Fl: AAM89250), Danio danio (zFF la: 
ААК54449), Danio danio (zfflb: AAF43283), Danio 
danio (zfflc: ААК19303), Danio danio (219: 
ААО59489), Drosophila | melanogaster (dFTZ-F1: 
P33244), Epinephelus  coioides (grouper ЕТ?-ЕІ: 
AAQ72771), Gallus gallus (CFTF/LRH-1: BAA22838), 
Gallus gallus (cSF-1: BAA22839), Homo sapiens (hSF-1: 
AAB53105), Homo sapiens (hFTF: AAD03155), 
Ictalurus punctatus (ccNR5AA: AAY45704), Mus 
musculus (Mlrh-1: AAA39447), Mus musculus (Msf-1: 
AAB28338), Oryzias latipes (mdFTZ-F1: BAA32394), 
Oncorhynchus | mykiss | (ttFTZ-Fl: ААМ83490), 
Oncorhynchus mykiss (rttLRH-1: BAE71417), Rana 
rugosa (rrFTZ-Fla: BAA94077), Rana rugosa (rrSF-1: 
BAA36789) and Trachemys scripta (tuSF-1: AAD01975). 


Real-time quantitative RT-PCR analysis of FTZ-F1 
expression in the tissues 

Total RNA was isolated from various tissues, СОМА 
was synthesized, and then real-time quantitative RT-PCR 
was conducted for analysis of FTZ-F1 expression using 
Roche lightcycler 480 Real-time PCR instrument (Roche, 
Switzerland) by means of the Real-time PCR Universal 
Reagent. A pair of gene specific primers (P5: 5'-ACACC 
CTTCCCACCATCTCTAA-3, P6: 5'-GGTGTAGTCGG 
GGCACTCG-3') were used to amplify a 194 bp FTZ-F1 
cDNA fragment by touchdown real-time PCR. A 211 bp 
p-actin fragment was amplified as an internal control 
with a pair of B-actin primers (P7: 5-GATGGTGGGTAT 
GGGTCAGAA-3', P8: 5-TTTGGGGTTCAGGGGAGC- 
3). Relative fluorescence unit, calculated threshold cycle 
(Ct), and dissociation curve were monitored by the 
analysis software of the system. Triplicate assays per 
RNA sample were carried out to determine the average 
Ct values. Data of relative expression levels were 
analyzed according to previously described methods 
(Livak & Schmittgen, 2001). All the PCR products were 
electrophoresed on 1.0% agarose gel and visualized 
using ethidium bromide staining. 


Statistical analysis 

The averages of the data obtained were subjected to 
one-way ANOVA and Duncan's new multiple-range test 
using SPSS 15.0 (Chicago, Ill, USA). Significant 
differences among treatment averages were determined 
ata P«0.05. 


RESULTS 


Isolation, characterization and phylogenetic analysis 
of the FTZ-F1 cDNA 

After RT-PCR and subsequent 5' and 3-RACE, a 1 
721bp FTZ-FIcDNA was obtained (Figure 1). The 
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СОС GGATCCACAGCCTACTG -260 
ATG ATCAGT СС АТОС AAA ACTCTG AGCTCTTG AATGGTC ATC ATCGT ACT. ACA CTTTCTA -240 
TTG GCC CTT CTGGCT CCTA A AGC ATCCCCGATAT AG AC AGCACTCACCCGG AATTC TOTI -180 
GGCCTT TCTTTC TTT ТТТ ATT TTG TTC ccc TCTGC AGCCCATGC AAGGCT GCTGCAACAC -120 
СТАСАТОСА GGTTTT AATC AT C AG ACAATCAT CACTG CCTGCC TGCTOTC T GCTGCT CAA -60 
ATGTTGGGAGACAAGGCTCACGGTGTGACTTTAAAGGTGATGGAGTACAC ATACGACGAG 60 


M L &G D К A H G V T E. HK cWo OMS GEO 95 СТ Y D E 20 


GACCTGGAAGAGCTGTGTCCTGTGTGCGGAGACAAAGTGTCTGGATATCA CTACGGCTTG 120 
D L E E L £t P wv (C G&G D K wv S RH Y MH Y G L 40 


CTCACCTGTGAGAGCTGCAAGGGCTTCTTCAAGAGGACAGTGCAGAACAA CAAGAGGTAC 180 
L T C E S € K G F F К R T М HG N N K R Y 60 


ACGTGTGCGGAAAACCAGGAGTGTATAATAGACAAGACCCAGAGGAAAAG ATGTCCCTTC 240 
T C A E N Q E С І | вк та в к R СОР F 80 


TGTCGGTTCCAGAAGTGCCTCAATGTCGGGATGCGGCTGGAAGCGGTGCG TGCAGACCGC 300 
t R F ü к с ІМ Уа м в Е A YSR A D R 100 


ATGCGTGGGGGTAGGAATAAATTTGGCCCCATGTACAAGAGAGACAGGGC CTTAAAACAG 360 
M в G G в м к ғ а P M Ү K R D В А L K Q 120 


CAGAAGAAAGCTTTGATACGATCCAACGGTTTCAAGTTGGAGAGCACGGC TCCTCCACCA 420 
а ққ A L | R S N G F K L E 8 T A P P P 140 


GCCTCTCCTCTGCAGGCAGACTACGGCTTCACTGGCACCTTGCACACCCT TCCCACCATC 480 
A 5 P L Q A D Y G F T G T сон T L P Т | 16 


TCTAAAAGCCTGCTCCCATCCACCCCGAGCTCCATCACTCCCACAGACTA CGAGGCTAAT 540 
Зх ЕСКЕК ЕЕ СЕН МКК ЛИН pur Be То SD Е A N 180 


CTCTACGGACCCCCGTCTTTGGGTATGGCCATGCAGTCACACGTGCCCCT CACCACCCAG 600 
L Y G P P S8 L а M À M Q S H ¥ P Е T T Q 200 


TACCAGTACACAGCCTTTCCCGGAAGGGCAATTAAAGCCGAGTGCCCCGA CTACACCAGC 660 
Y Q Y T A F P G В А 1 К А E C P D *Y T $ 220 


TCCCCTGAGACTCTTACAGGATACCCTTACCCAGACATGTACCCCTCAGC CTCGCCGCAA 720 
Se GS СЕТ, L T G Y P Y P D M Y P 5 A S P о 240 


CCCCCCAGCCTGCCGCCGCTGGTTCTGGAGCTGCTGCGCTGTGACCCGGA CGAGCTGGAC 780 
РР Sr OL GR Pe Cb. Wo LE? КЖ © бов E L D 260 


GTCCAGAACAAGATTGTGACTCACCTGCAACAGGAGCAAAACGGCCGGGG TCGGCTGGAA 840 
мо N K І У T H L QQA E Q N G R G R L E 280 


AAGCCCAGCACCTTCAGCCTCATGTGTCGCATGGCAGATCAAACGTTGTTCTCCATCGTG 900 
Ke PAST FP Se d. SMe € ROOM. А CD. GT. LAFE $ м 300 


GAGTGGGCTCGGAGCTGCATCTTCTTTAAGGAGCTGAGGGTGGGAGATCA AATGAAGCTS 960 
E WoA R $ E | Е-Е K E L R МБЕ D M K L 320 


CTTCATAACTGCTGGTCTGAACTGCTGGTCCTGGATCACATTTTCAGACA AGTGCAACAC 1020 
L H N € W S EL L v L DH | F ВО ¥ G H 340 


GCAAAGGAAGACAGCATTCTGCTGGTGACCGGCCAGGAGGTGGAGCTGTC ATCCATCCTG 1080 
A K E D S | LLY T G6 Q E V E L $ | L 360 


TCTCAGGCTGAGGGGACACTCTCCAGTCTGGTCCAAAGAGGTCAGGAGCT GGCAGCGAGG 1140 
8 Q 4 E G T L 8 38 L V @ R G Q E L ^ A 350 


CTGCGGGTGCTGCAGGTCGACCGCAGAGAGATCGCCTGTCTGAAATTCCTCCTCCTTTTC 1200 
L R ¥ L Q ¥ D R R E | А С L K Е L E, ES CHE 400 


AACCCAAATGTAAAATTGCTGGAGAACCAGGCGTTTGTGGAGGGCGTCCA GGAGCAGGTG 1260 
МОР WN УК L LE НОА F У E Б ¥ @ Е Я Vv 420 


AATGCTGCTCTGCTGGAGTACACCCTGTCTGCCTACCCTCAGTTCCAAGA GAAGTTCAGC 1320 
N A A L L E Y T L S8 A Y P QO F Q E КЕ S 440 


CAGCTGTTGGTACGGTTGCCGGAGCTGCGCTCCCTCAGCACGCAGGCCGA GGACTACCTA 1380 
Qo oL oL oN RD ЕР Б 5 15027 20: OE D Y L 460 


TGCTACATGCACGTGAGTGGAGAGGTGCCCTGCAACAACCTGCTGATCGA GATGCTGCAC 1440 


U Y M H ¥ S G E УР C N N L L | E M L H 480 
1461 
GCTAAGAGGGCGTGTGTGTIGA 486 


A K R A € V 


Figure 1 cDNA nucleotide and predicted amino-acid sequences of the GIFT Oreochromis niloticus FTZ-F1 
Nucleotide are indicated above and numbered to the right of each lane (upper row). The deduced amino acid sequence is shown below the nucleotide sequence. 


Amino acids are number to the right of each lane (lower row). 


FTZ-F1 cDNA included an open reading frame, bind regions (LII and III, FTZ-F1 box) and the activation 
encoding a predicted 486 amino acid protein, and function-2 hexamer (Figure 2).The present sequence is 
including the highly conserved DNA-binding and ligand- available in GenBank database with the accession 
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giFTZ-Fi ——--MLGDEAHGVTLEVMEYTYDEDLEELCPYCGDEVSGYHYGL 40 
zdFTZ-Fi -MLGDKAHGVTLEVMEYTYDEDLEELCPVCGDEVSGYHYGL 40 
ЖЕНІП) 0—————— — —-—-MLGEKAQGVPLEVMDYTYDADLEELCPVCGDEVSGYHYGL 40 
сЕТЕ/І.ЕН-1 MLPEVETEALGLARSNGEQGQMPENMQVSQFEMVNY SYDEDLEELCPYCGDKVSGYHYGL 60 
rrFIZ-Fia MLPEVESEALGLSRSHGEHGQMPDNMQVSQFKMMCYSYDDDLEEMCPVCGDKVSGYHYGL 60 
zFFla MLPEVESEYLGLARSHGEQGHMPGNMQAPQFKEMMDYSYDEDLDEMCPVCGDEVSGYHYGL 60 
hSF-1 ———--MDYSYDEDLDELCPYCGDEVSGYHYGL 27 
Msf-1 ——---MDYSYDEDLDELCPVCGDKVSGYHYGL 27 
grouper ———-— — —MEHRHDVGLEELCFVCGDEVSGYHYGL 27 
1211Ж KDK Doloiokdolololookolcialcial 
giFIZ-Fi LTCESCKGFFERTVQNNERYTCAENQECIIDKTQRERCPFCRFOKCLNVGMRLEAVRADR 100 
zdFTZ-Fi LTCESCKGFFERTVONNERYTCAENQECEIDRTQRKRCPFCRFOECLNVGMRLEAVRADR 100 
rtiFTZ-Fi LTCESCKGFFERTVQNNERYTCAENQECKIDKTQRERCFFCRFQKCLNVGMRLEAVRADR 100 
cFTF/LRH-1 LTCESCKGFFKRTVQNNKRYTCIENQNCQIDKTQRKRCPYCRFQKCLSVGMKIEAVRADR 120 
rrFIZ-Fia LICESCKGFFERTVOQNNERY TCI ENQTCQIDETQRERCPYCRFQECLSVGMELEAVRADR 120 
zFFla LTCESCEGFFERTVQNNERYTCIENQSCOIDKTQORKRCPYCRFOKCLTVGMELEAVRADR 120 
hSF-i LICESCKGFFERTVQNNEHY TCTESQSCE IDETQRERCPFCRFQECLTVGMRLEAVRADR 57 
Msf-1 LICESCKGFFERTVONNEHY TCTESQSCE IDETQRERCPFCRFQECLIVGMRLEAVRADR 87 
grouper LTCESCKGFFERTVQNNKEY ICAEKQDCRIDITQRKRCPFCRFQKCLHVGMRLEAVRADR 87 
ЖЖЖ oook o Ж Ж. Ж ж ЖЖ НОЕ С ЖЖЖ жж жж 
giFTZ-Fi MRGGRNEPGPMYERDRALEOQQKEAL IRSNGFELESTAPPPASPLQADYGFTGILHTLFTI 180 
zdFIZ-Fi MRGGRNEFGPMYERDEALEDQEKAL IRSNGFELESTALPSASPLQTDYSFTGILHTLFTI 180 
rtFTZ-Fi MRGGRNKFGPMYERDRAMERGEKALIRSTGFKLDS-APPQISPVQTNYGFTGTLHSLPSL 159 
cFTF/LRH-1 MRGGRNEFGPMYERDRALEGQOEKAL IRANGLKLEA-MTQVIQAMPTDLTISSAIQNIHSA 179 
rrFIZ-Fia MRGGRNEFGPMYERDRALEDQEKALIRANGLKLEA-MSQVIQAIPTDLITISSAIQNIHSA 179 
zFFia MRGGRNKFGPMYERDRALEDQEKALIRANGLKLEA-MTQVMQTVPADLTITSAIQNIHSA 179 
hSF-1 MRGGRNKFGPMYERDRAI EGGQERAQTRANGFKLET-GPPMGVPPPPPPAPDYVLPPSLHG 146 
Msf-1 MRGGRNKFGPMYERDRALEGQKEAQIRANGFELET-GPPMGVPPPPPPPPDYMLPPSLHA 146 
grouper MRGGRNEFGETYELDEALEI QREAL IQASGFRLDSNPALV: SSTHQRDLT FTGDLHLAPIL 147 
жасаа ЖЖ dole oka RD DKS СЖ: 
giFIZ-Fi SESLLPSTP—— ———-—— SSITPIDYEANLYGPPSLGMAMQSHV-PLTTQYQYTA 205 
zdFTZ-Fi SESLLPSTT— SSVIPTDYEANLYGPTSLGMAMQSHV-PLNPQYQYTA 205 
rtFIZ-Fi PEGLMPPIP— TSINPIDYEASLYGPPSLGVAMQPNG-PLPTQYQYTA 204 
cFIF/LRH-1 SKGLPLNHT— ALPPTDYDRSPFVISPISMIMPPHG-SLOGYQTYGH 223 
rrFIZ-Fia -ALPPTDYDRSPFVTSPISMAMPPHG-SLQSYQAYGH 223 





zFFla SRGLPLSHEHHEHHHHEHIHSSSS АІ PPADFDRSPFVTSPY SMAMPPHAGGLQGYQAYGH 239 

















hSF-i —PPAGP-LGDFGAPALPMAVPGAHGPLAGY-LYPA 187 
Msf-1 —PPSGP-LGDFGAPSLPMAVPGPHGPLAGY-LYPA 187 
grouper ——LPSTQNSHMSYQPPSLCSLLPSSS-PGATQYQCIS 185 
gfFIZ-Fi FPGRAIKAECPD-YTSSPE-TLIGYPYPDMYPS-ASPQPPSLPPLVLELLRCDPDELDVQ 262 
zdFIZ-Fi FPSRAIKAECPD-YISSPE-SLTGYPYPDVYPS-ASPQPPSLPPLVLELLRCDPDELVWQ 262 
riFIZ-Fi FPSRAIKSEYPDPYTSSPE-SLVGYPLPDGYPSGGSPQLPSQPPLVLELLRCDPDEMQVQ 263 
cFTF/LRH-1 FPSRAIKSEYPDPYTSSPE-SIMGYSYMDGYQT-—-SSPASIPHLILELQKCEPDEPQVQ 279 
rrFIZ-Fia FPSRAIKSEYPDPYTSSPE-SLMGYSYMDSYQS—-SSPSNIPHLIVELLKCEPDEPQVQ 279 
zFFia FQSRTIKSEYPDPYTSSPE-SLMGYPYVEAYAG-——-GSPPSFPHLVYELLKCEFDEPQVQ 295 
hSF-i FPGRAIKSEYPEPYASPPQ-PGLPYGYPEPFSG-—G—-FNVPELILQLLQLEPDEDQVR 241 
Msf-1 FSNRTIKSEYPEPYASPPQQPGPPTSYPEPFSG— G—FNVPELILQLLQLEPEEDQVR 242 
grouper ee ETAG NSOE re eee eee 244 
- ІжІЖ тж. ж а 2. Ж > DOIXIX* 
giFIZ-Fi NKIVIHLQQEQNGRGRLEKPSTFSLMCRMADGTLFSIVEWARSCIFFKE]R 
zdFTZ-Fi NKIVAHLQQEQSNRGRLDEPSTF SLMCRMADQTLFSIVEWARSCIFFKELR 
riFIZ-Fi SKIMAFLOQQEQSGRGRQEEL STF SLMC RMADQTLFSIVEWARSCIFFKELE 
сЕТЕ/ІЕН-і AKIMAYLQQEQANRSKHEKLNTFGLMCKMADQILF5 [VEWARSSIFFREPR 
rrFIZ-Fia SKIMSYLQQEQANRSKHDELNTFGLMC KMADGTLFS IVEWARSSIFFRELE 
zFFia ARILAYLQQEQASRGKHEKLNTFGLMC EMADQTLFSIVEWARSSIFFREER 
hSF-i ARILGCLQEPTKSRP--DQPAAFGLLCR 
Msf-1 ARIVGCLOEPARSRS реРАРРЭ САТТИ 
grouper TQITARLLQEHMSWEEHGNP STF SLMC. K 
IO e n E І ж ЖЖ XXORX 7 DXODRIEODK 2 
gfFTZ-F1 NCWSELLVLDHIFRQVQHAKEDSILLVIGQEVELSSILSQAEGTLSSLVQRGQELAARLR 382 
zdFIZ-Fi NCWSELLVLDHIFRQVQHGREDS ILLVTGQEV@LSSILSQGEATLCSLVQRGQELAARLR 382 
xtFTZ-Fl NCWSELLVLDHVFRQVQHGEESSLLLVIGQEMDLSSMGSQAGVILSGLVQRGQVLAGRLL 383 
cFTF /LRH-1 NCWSELLILDHIYRQVVHVKEGSILLVTGQQVDYSVIASQAGATLNNLMSHAQELVAKLR 399 
xrFIZ-Fia NCWSELLILDHIFRQVLHGKEGSILLVTGQQVDFSVIVSQAGTTLNNLMSHAQELVARLR 399 
zFFla NCWSELL ILDHVFRQVMHAKEGS ILLVIGQQVDYALIASQAGAT LNNLLSHAQELVSELR 415 
hSF-i NCWSELLVF DHIYRQVQHGKEGS [LLVIGQEVELTIVATQAGSLLHSLVLRAQELVLQLL 359 
Msf-1 NCWSELLVLDHIYRQVQYGKEDSILLVIGQEVELSTVAVQAGSLLHSLVIRAQELVLQLH 360 
grouper SCWSELLLLDIISRQVLCGKEGSLLPFTGQEMDLSDMASHAGLTLASLVQRGQELVEKLH 364 
.OOEWOROR I OK 1 ЖЖЖ ЖЖ. kik Jaa 1 І = ж .x* :1.ж ж. ІЖ 
g£FIZ-Fi VLOVDRREIACLEFLLLENPN/KLLENQAFVEGVQEQVNAALLEYTLSAYPQFQEKFSQL 442 
zdFIZ-Fi ALOVDREEIACLEFLLLENPN/ELLENQAFVEGVQEQVNGALLEYTLTTYPQFQERFSQL 442 
xtFTZ-Fl TLQVDRREVACLEFLLLFNPNVELLENQAFVESVQEQVNGALLEYTLCTYPLYLDEFSQV 443 
cFTF/LRH-1 SLQFDLREFVCLEFLVLF SLINKNLENFQLVEGVQEQVNAALLDYTMCNYPQQTDKFGQL 459 
rrFIZ-Fia SLGFIMREFVCLEFLVLFSLIVENLENFQLVEGVQEQINAALLDYTLCYYPQQTDKFGQL 459 
zFFla SLOLDORFFVCLEFLVLFSLINENLENFHLVESVQEQVNAALLDYVMCNYPQQTDKFGQL 475 
hSF-i ALQLDROEFWCLEFITLFSLD KFLNNHILVKDAQEKANAALLDYTLCHYPHCGDKFQQL 419 
Msf-1 ALQLDRQEFWCLEFLILFSLINEFLNNHSLVKDAQEKANAALIDYTLCHYPHCGDKFQQL 420 
grouper ILEVURGEFACIKFLILENPDNTPLEDHQLIESVOEQVEGALLEYTLCTSSHPPGHFAHL 424 
LL DX LRRD eR 11. ЖОС III.. kI 1. Жж Ж nt 
gfFTZ-Fi LVRLPELRSLSTQAEDYLCYMHVSGEVPCNN LLIEM НАКРАСУ 486 
zdFIZ-Fi VVRLPELRSLSTQAEDYLCYMHL SGEVPCNN // ГЕМ НАКЕАСУ 486 
rtFIZ-Fi VMRLPELRALSTQAEDYLCTKHLSGEVPCNN LLIEM HAKRACV 487 
cFTF/LRH-1 LLRIPEIRAISMQAEEYLYCKHLNGDVPCNNILIEM НАКВА— 501 
rrFIZ-Fia LLRLPEIRAISLQAEEYLYYKHLNGDVPCNN LrIEM-HAKRA— 501 
zFFia LLRLPEIRAISLQAEEYLYYKHLNGDVFCNN LLIEM-HAKRA— 517 
hSF-i LLCLVEVRALSMQAKEYL YHRHLGNEMPRNN LITEM QARQT— 461 
Msf-1 LLCLVEVRALSMQAKEYLYHKHLGNEMPRNNLLIEM QAKQT— 462 
grouper LLCLSELRCLSSFAEDYLYCKHLSGEAPCNN LrTEM-HAKHSWA 468 
To Жж жж Ж:..: Ж ЖЖЖ fick Ek. 


Figure 2 Amino acid alignment of the GIFT tilapia FTZ-F1(gfFTZ-F1) with other species 
Oryzias latipes (mdFTZ-F 1: BAA32394), Oncorhynchus mykiss (rtFTZ-F1: AAW83490), Gallus gallus (CFTF/LRH-1: BAA22838), Rana rugosa (rrFTZ-Fla: 
BAA94077), Danio danio (zFF 1a: AAK54449), Homo sapiens (hSF-1: AAB53105), Mus musculus (Msf-1: AAB28338), Epinephelus coioides (grouper FTZ- 
ЕІ: AAQ72771). The highly conserved regions of I ,II,III, FTZ-F1 box and AF-2 hexamer are shown іп gray fonts, blace letters in black boxes, gray and 
boldface fonts, boldface fonts, and italics, respectively. 
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number JF742993. 

The deduced amino acid sequence of the GIFT 
tilapia FTZ-F1 had 39%-94% of identity with that of 
other vertebrates. A higher level of identity was found 
when the GIFT tilapia FTZ-F1 was compared to the 
other teleosts FTZ-F1 proteins. The deduced amino acid 
sequence of the GIFT tilapia FTZ-F1 had 94%, 93%, 
83%, 81% identity with that of the O. latipes FTZ-F1, A. 
schlegelii FTZ-Flb, О. mykiss FTZ-F1, С. semilaevis 
FTZ-F1, respectively; 68%, 68%, 65%, with that of the Т 
scripta SF-1, G gallus 8Е-1, В. rugosa SF-1, 
respectively; 63%, 62% with that of G gallus LRH-1, M. 
musculus LRH-1, respectively. The lowest sequence 
identity was with the D. melanogaster FTZ-Fl. The 
GIFT tilapia FTZ-Fl was 68%, 67%, 63%, 46% 
sequence identity with D. danio fflb, ffld, ffla, fflc, 
respectively. The GIFT tilapia FTZ-Fl amino acid 
sequence and several published FTZ-F1 (as shown іп 
Figure 3) in various vertebrate species were used to infer 
phylogenetic relationships. The phylogenetic tree shows 
that the GIFT tilapia clustered firmly with other teleosts 
FTZ-F1 (Figure 3). 


Tissues expression of FTZ-F1 

The ЕТ?-ҒІ mRNA expression in different tissues 
of the GIFT tilapia was analyzed using real-time PCR. 
The FTZ-F1 transcripts were highly expressed in testis 
and ovary and weakly in kidneys (Figure 4). No 
transcripts were found in brain, heart, liver, spleen, 
intestine, gills and muscle. The amounts of FTZ-F1 
mRNA expression in the testis of adult male fish were 
approximately 3 fold of that in the male kidney (Figure 4) 
and the amounts of FTZ-Fl mRNA expression in the 
ovary of adult female fish were approximately 2.78 fold 
of that in the female kidney (Figure 4). 


DISCUSSION 


We successfully isolated a homologues FTZ-F1 
from the GIFT tilapia ovary. Amino acid alignment 
analysis showed that the FTZ-F1 was characteristic of 
nuclear receptor superfamily, with highly conserved 
regions of Г,П,Ш, ЕТ?-ЕІ box and AF-2 hexamer. 
Phylogenetic analysis of the GIFT tilapia FTZ-F1 
showed that it clustered with other teleost FTZ-F1. 

Likewise, the GIFT tilapia FTZ-F1 was expressed in 
the gonad and kidney, which suggested that the FTZ-F1 
belonged to SF-1/Ad4BP (Hammer & Ihgraham, 1999; 
Parker & Schimmer, 1997). No transcripts of the GIFT 
tilapia FTZ-F1 were found in the liver, which implied 
that the FTZ-F1 did not belong to the LRH/FTF group 
(Lu et al, 2001). 

There are a variety of evidences for FTZ-F1 activity 
in other fish. During embryogenesis of zebrafish, the 
expression of homologues FTZ-F1 was detected in the 
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pituitary, mandibular arch, pronephric duct, liver, rostral 
diencephalons, hindbrain, and pancreas (Chai & Chan, 
2000; Von Hofsten et al, 2001), indicating that they 
might be involved in tissue differentiation. The tongue 
sole FTZ-F1 mRNA was detected in all stages from 
zygote to 25 days after hatching of the tongue sole, 
suggesting that FTZ-F1 might be involved in the 
organogenesis of the tongue sole embryo (Deng et al, 
2008). The orange-spotted grouper E. coioides is a 
protogynous hermaphroditic fish, and the expression of 
FTZ-F1 in the gonad also decreased significantly in 
response to MT treatment (Zhang et al, 2007). 
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Figure 3 Neighbor joining tree of the GITF Tilapia FTZ-F1 
апа FTZ-F1 from other organisms based on amino 


acid sequences 

Distances are used to construct the phylogenetic tree and bootstrap values 
based оп 1 000 resampling replicates. The bottom scale refers to percentage 
divergence (p-distance). Acanthopagrus schlegelii (55ЕТ7-Е1а: AAS75791), 
Acanthopagrus schlegelii (s)FTZ-F 1b: AAS75792), Bos taurus (btAd4BP: 
BAA02764), Cynoglossus semilaevis (tsFTZ-Fl: ABQ41307), Clarias 
gariepinus (acFTZ-Fla: AAG49004), Clarias gariepinus (acFTZ-F1b: 
AAG49005), Carassius auratus (Gftz-Fl: AAM89250), Danio danio 
(zFFla: AAK54449), Danio danio (zfflb: AAF43283), Danio danio (zfflc: 
AAK19303), Danio danio (zffld: AAO59489), Drosophila melanogaster 
(dFTZ-F1: P33244), Epinephelus coioides (grouper FTZ-F1: AAQ72771), 
Gallus gallus (cFTF/LRH-1: BAA22838), Gallus gallus (cSF-1: 
BAA22839), Homo sapiens (hSF-1: AAB53105), Homo sapiens (hFTF: 
AAD03155), Ictalurus punctatus (ccNR5AA: AAY45704), Mus musculus 
(Mirh-1: AAA39447), Mus musculus (Msf-1: AAB28338), Oryzias latipes 
(mdFTZ-F1: BAA32394), Oncorhynchus mykiss (rtFTZ-F1: AAW83490), 
Oncorhynchus mykiss (rtLRH-1: BAE71417), Rana rugosa (rrFTZ-Fla: 
BAA94077), Rana rugosa (rrSF-1: BAA36789), GIFT tilapia CgfFTZ-Fl: 
1Е742993) and Trachemys scripta (tuSF-1: AAD01975). 
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Figure 4 FTZ-F1 expression from various tissues of adult 
female and male, GIFT tilapia 
1: testis; 2. brain; 3. heart; 4. kidney; 5. spleen; 6. intestine; 7. 
Sgills; 8. muscle; 9. Qovary; 10. Qbrain; 11. Qmuscle; 12. Qheart; 13. 
Oliver; 14. Qkidney; 15. Qspleen; 16. intestine; 17. Qgills; 18. liver. 
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Preliminary analysis of the mitochondrial genome 
evolutionary pattern in primates 
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Abstract: Since the birth of molecular evolutionary analysis, primates have been a central focus of study and mitochondrial DNA is 
well suited to these endeavors because of its unique features. Surprisingly, to date no comprehensive evaluation of the nucleotide 
substitution patterns has been conducted on the mitochondrial genome of primates. Here, we analyzed the evolutionary patterns and 
evaluated selection and recombination in the mitochondrial genomes of 44 Primates species downloaded from GenBank. The results 
revealed that a strong rate heterogeneity occurred among sites and genes in all comparisons. Likewise, an obvious decline in primate 
nucleotide diversity was noted in the subunit rRNAs and tRNAs as compared to the protein-coding genes. Within 13 protein-coding 
genes, the pattern of nonsynonymous divergence was similar to that of overall nucleotide divergence, while synonymous changes 
differed only for individual genes, indicating that the rate heterogeneity may result from the rate of change at nonsynonymous sites. 
Codon usage analysis revealed that there was intermediate codon usage bias in primate protein-coding genes, and supported the idea 
that GC mutation pressure might determine codon usage and that positive selection is not the driving force for the codon usage bias. 
Neutrality tests using site-specific positive selection from a Bayesian framework indicated no sites were under positive selection for 
any gene, consistent with near neutrality. Recombination tests based on the pairwise homoplasy test statistic supported complete 
linkage even for much older divergent primate species. Thus, with the exception of rate heterogeneity among mitochondrial genes, 
evaluating the validity assumed complete linkage and selective neutrality in primates prior to phylogenetic or phylogeographic 
analysis seems unnecessary. 


Keywords: Mitochondrial genome; Evolutionary pattern; Codon usage bias; Complete linkage; Evolution neutrality; Primates 


The mitochondrial genome of vertebrates exhibits 
several peculiar features including maternal inheritance, 
the presence of single-copy orthologous genes, a lack of 
recombination, evolutionary neutrality, and a high 
mutation rate—characteristics that make it well suited for 
evolutionary studies (Pesole et al, 1999; Saccone et al, 
2000). Among all the characteristics and assumptions of 
the mitochondrial genome, the two most important 
aspects for explaining the evolutionary process are the 
presumed freedom from positive Darwinian (adaptive) 
natural selection and lack of recombination. 
Mitochondrial proteins are central to the cellular 
oxidative phosphorylation pathway and are functionally 
conserved across metazoan phyla (Gray et al, 1999). The 
importance of these mtDNA protein products supports 
the hypothesis that some variation in the mtDNA 
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molecule is adaptive. Nonetheless, of the numerous 
single mtDNA gene studies and several complete 
mitochondrial genome comparisons published over the 
past two decades, very few have been able to reject 
neutrality in favor of positive Darwinian selection 
(Garcia-Martinez et al, 1998; Rand et al, 1994). Indeed, 
previous studies on vertebrates indicated that different 
levels of variability are attributable to functional 
constraints and/or slightly deleterious polymorphisms, a 
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pattern consistent with near neutrality (Nachman et al, 
1996; Templeton, 1996). 

Recombination breaks down the correlation in 
genealogical history between different regions of a 
genome, which may be the incorrect inference of 
evolutionary history (Schierup & Hein, 2000). The issue 
of whether recombination occurs in the human 
mitochondrial genome remains controversial as the 
necessary enzymes for recombination are in fact present 
in the mitochondria, and a few paternal mitochondria do 
penetrate the egg during fertilization (Thyagarajan et al, 
1996), providing evidence that suggests recombination is 
possible, at least in humans. Recent broad surveys of 
animal mitochondrial genomes have also concluded that 
recombination is widespread (Piganeau et al, 2004; 
Tsaousis et al, 2005). Another interesting issue about 
evolutionary patterns of the mitochondrial genome is the 
codon usage bias in its protein-coding genes. Many 
studies have analyzed the differences in codon usage 
between different organisms and between proteins within 
the same organism (Bulmer, 1991; Kanaya et al, 2001; 
Sharp et al, 1995) and one of their main conclusions is 
that there are strong correlations between codon usage 
and genomic GC content. Moreover, human codon usage 
may be determined solely by GC content апа its 
isochores composition (Kanaya et al, 2001). Despite 
suggestive evidence and hypothesis, the relationships 
between codon usage and the underlying evolutionary 
constraints are still not fully understood (Prat et al, 2009). 

Since the birth of molecular evolutionary analysis, 
the evolutionary relationships of our own order, primates, 
have been of central interest. Strangely, no comprehensive 
and accurate evaluation for the nucleotide substitution 
pattern(s) of the various mtDNA genes (1.е., tRNA and 
rRNA genes, synonymous and nonsynonymous positions 
of protein-coding genes) have been conducted to date. 
Most primate studies demonstrated mtDNA evolutionary 
neutrality and no recombination, and generally reported 
mtDNA still evolves about 5- to 10-fold more rapidly 
than single-copy nuclear DNA (nDNA) (Brown et al, 
1982) at an overall nucleotide substitution rate of 1% per 
million years (Brown et al, 1979; Wilson et al, 1985). 

Аз large libraries of whole-mtDNA genome 
sequences from primates accumulate, we are able to 
conduct more sophisticated genome-wide analyses to 
obtain solid information on the pattern of whole-mtDNA 
genome evolution in primates while also determining the 
relative evolutionary rate of each mitochondrial 
component. Here, we report on the rates and patterns of 
mtDNA sequence variation and the codon usage in 
primates and describe the results of tests conducted to 
detect recombination and selection. These results will not 
only be of interest in their own right but will also have 
implications for the use of mtDNA in evolutionary 
studies of primates. 
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Materials and Methods 


Sequences 

54 Mitogenomic sequences from 44 primate species 
were downloaded from GenBank (Table 1). All 54 
sequences were aligned, edited, and compared using the 
Sequencher 4.5 (Gene Codes Corporation, Ann Arbor, 
MI). We excluded all gaps and ambiguous alignment 
sites, resulting in 14,764 bp of sequence from the 13 
protein-coding, 2 rRNA, and 21 of the 22 tRNA genes. 
Protein-coding sequences from the L-strand-encoded 
genes (ND6 and 8 tRNA genes) were converted into 
complementary strand sequences. The tRNA-Glu gene 
and the Control Region were not included because the 
alignments of those sequences are difficult due to a large 
number of highly variable sites, insertions, and deletions. 


Phylogenetic analyses and sequence diversity 

The genealogical relationships between the 
mitochondrial genomes were analyzed by parsimony 
using PAUP* (Swofford, 2002). Bootstrapping (Felsenstein, 
1985) was used to test monophyly. For this study, 1,000 
pseudosamples were generated to estimate the bootstrap 
proportions. 

Nucleotide diversity (п) was calculated for all 13 
protein-coding, 2 rRNA, and 21 concatenated tRNA 
genes. Measures of sequence diversity were also 
calculated by the sliding window method using the 
DnaSP 4.10 (Rozas et al, 2003). A window (500 bp in 
length) was moved along the sequences in steps of 50 bp. 
Patterns of nucleotide diversity for all positions as well 
as patterns of synonymous and nonsynonymous 
nucleotide variation were calculated in each window, and 
the value was assigned to the nucleotide at the midpoint 
of the window. The gamma parameter alpha (а), which 
represents the extent of rate heterogeneity among sites, 
was estimated from the individual data sets with TREE- 
PUZZLE 5.2 (Schmidt et al, 2002) using the 8-categories 
option. 


Codon usage bias 

Several parameters related to codon usage bias, 
such as the codon bias index (Morton, 1993), the 
effective number of codons (Wright, 1990), and G + C 
content at second and third positions as well as overall 
were estimated for each primate mitochondrial protein- 
coding genes using DnaSP 4.10 (Rozas et al, 2003). 
Furthermore, to determine whether the compositional 
changes of the nucleotide content in the primate 
mitochondrial protein-coding genes are caused by 
directional mutation pressure or a result of positive 
selection, we performed a correlation analysis. If the 
nucleotide bias affects both the synonymous and 
nonsynonymous sites in protein-coding genes, positive 
selection could not solely explain this phenomenon, 
because positive selection should not affect silent 
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Table 1 Primate taxa and GenBank accession numbers for mitochondria genome genes. 
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Galagidae Galago senegalensis AB371092 Presbytis melalophos DQ355299 
Otolemur crassicaudatus AB371093 Pygathrix nemaeus DQ355302 
Lorisidae Nycticebus coucang AJ309867 Pygathrix roxellana DQ355300 
Loris tardigradus AB371094 Nasalis larvatus DQ355298 
Perodicticus potto AB371095 Chlorocebus pygerythrus 1 EF597501 
Daubentoniidae Daubentonia madagascariensis 1 AM905039 Chlorocebus pygerythrus 2 EF597500 
Daubentonia madagascariensis 2 AB371085 Cercopithecus aethiops sabaeus DQ069713 
Indriidae Propithecus verreauxi AB286049 Papio hamadryas Y18001 
Lemuridae Varecia variegata AB371089 Theropithecus gelada FJ785426 
Lemur catta AJ421451 Macaca sylvanus AJ309865 
Eulemur macaco AB371088 Macaca thibetana EU294187 
Eulemur mongoz AM905040 Macaca mulatta AY612638 
Eulemur fulvus mayottensis AB371087 Macaca fascicularis FJ906803 
Eulemur fulvus fulvus AB371086 Hylobatidae Hylobates lar X99256 
Tarsiidae Tarsius syrichta AB371090 Pongo pygmaeus D38115 
Tarsius bancanus AF348159 Pongo abelii X97707 
Cebidae Saguinus oedipus FJ785424 Gorilla gorilla 1 D38114 
Cebus albifrons AJ309866 Gorilla gorilla 2 X93347 
Aotidae Aotus lemurinus FJ785421 Pan troglodytes 1 EU095335 
Atelidae Ateles belzebuth FJ785422 Pan troglodytes 2 D38113 
Pitheciidae Callicebus donacophilus FJ785423 Pan troglodytes 3 X93335 
Cebidae Saimiri sciureus 1 AB371091 Pan paniscus D38116 
Saimiri sciureus 2 FJ785425 Homo sapiens 1 АМ948965 
Cercopithecidae Procolobus badius DQ355301 Homo sapiens 2 Х93334 
Colobus guereza AY863427 Homo sapiens 3 D38112 
Semnopithecus entellus DQ355297 Homo sapiens 4 AM711903 
Trachypithecus obscurus AY863425 Homo sapiens 5 AF346992 


nucleotide positions. The correlation analysis of GC 
content at the second codon position (nonsynonymous 
mutations related, GC?) and GC content at third position 
(synonymous mutations related, СС;з) with CG content 
of all protein codon genes (GCc) were implemented 
using R 2.6.2 (Ihaka & Gentleman, 1996). 


Tests of neutrality 

To investigate protein evolution, we first calculated 
dy/ds ratios among the set of 54 primate data and tested 
their significance using Z-tests (Nei & Kumar, 2000) for 
each of the 13 protein-coding genes as implemented by 
MEGA 4.0 (Tamura et al, 2007). We also tested for site- 
specific positive selection with Bayesian posterior 
probabilities under the Ny98 fitness regime (0<@, <1, 
05-1, and оз 71) (Nielsen & Yang, 1998) as described 
and implemented by MrBayes 3.1.2 (Huelsenbeck & 
Ronquist, 2001; Ronquist & Huelsenbeck, 2003). 


Tests of recombination 
We calculated the Kelly (1997) measure of linkage 
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disequilibrium (ZnS); ZnS being the average of В? (Hill 
& Robertson, 1968) over all pairwise comparisons. We 
also employed the pairwise homoplasy test (Фу ), а 
relatively new and robust statistical estimator of 
recombination (Bruen et al, 2006) that is a coalescent- 
based estimator of the genealogical correlation or 
compatibility among sites, which is negatively correlated 
to recombination. The estimator was calculated, and its 
significance was estimated via a permutation test using 
Splitstree 4.10 (Huson & Bryant, 2006). 


Results 


Phylogenetic analyses and sequence diversity 

The parsimony tree for primates, based on the 
mitochondrial genome, is shown in Figure 1. Each 
internal node was supported by more than 80% bootstrap 
(data not shown). The mtDNA polymorphism data and 
estimates of nucleotide divergence are summarized in 
Table 2, and the distributions of sequence diversity 
across the mitochondrial genome are presented in Figure 
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Galago senegalensis 
Otolemur crassicaudatus 
Nycticebus coucang 
Loris tardigradus 
Perodicticus potto 
Propithecus verreauxi 
Varecia variegata 

Lemur catta 

Eulemur macaco 
Eulemur mongoz 
Eulemur fulvus mayottensis 
Eulemur fulvus fulvus 
Daubentonia madagascariensis 1 
Daubentonia madagascariensis 2 
Tarsius syrichta 

Tarsius bancanus 
Saguinus oedipus 

Aotus lemurinus 

Cebus albifrons 

Saimiri sciureus 1 
Saimiri sciureus 2 

Ateles belzebuth 
Callicebus donacophilus 
Procolobus badius 
Colobus guereza 
Semnopithecus entellus 
Trachypithecus obscurus 
Presbytis melalophos 
Nasalis larvatus 
Pygathrix nemaeus 
Pygathrix roxellana 
Papio hamadryas 
Theropithecus gelada 
Macaca sylvanus 

Macaca thibetana 
Macaca mulatta 

Macaca fascicularis 
Chlorocebus pygerythrus 1 
Chlorocebus tantalus 
Chlorocebus pygerythrus 2 
Hylobates lar 

Pongo abelii 

Pongo pygmaeus 

Gorilla gorilla 2 

Gorilla gorilla 1 

Pan paniscus 

Pan troglodytes 1 

Pan troglodytes 2 

Pan troglodytes 3 

Homo sapiens 1 

Homo sapiens 4 

Homo sapiens 5 

Homo sapiens 3 

Homo sapiens 2 


Figure 1 Parsimony tree of 54 primates based on mitochondrial genome (H;- Н, were selected nodes for recombination tests) 


2. An obvious decline in primate nucleotide diversity was 
noted in the subunit rRNAs and tRNAs as compared to 
the protein-coding genes. Among the 13 protein-coding 
genes, ATP8-ATP6 and the NADH dehydrogenase (ND) 
complex (with the exception of ND1) showed higher 
rates of sequence diversity, while the cytochrome 
oxidase (COX) complex exhibited lower diversity within 
Primates. Notably, the nonsynonymous divergence plots 
were similar to those obtained from the total nucleotide 
divergence, but the trend was not observed in the sliding- 
window plot of synonymous divergence. Although the 
COX complex exhibited an obvious decline in total 
nucleotide divergence and nonsynonymous divergence, 
these genes had higher rates of synonymous diversity. 
Based on the estimated gamma parameter alpha (o), the 
rate heterogeneity among sites was not even, but the 
pattern of heterogeneity among sites was distributed 
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across the primate mitochondrial genome. 


Codon usage 

Several parameters related to codon usage bias were 
estimated to check whether synonymous mutations are 
selectively neutral. The codon bias index (CBI) is a 
measure of the deviation from the equal use of 
synonymous codons with values ranging from 0 (uniform 
use of synonymous codons) to 1 (maximum codon bias). 
CBI values were intermediate in our study, ranging from 
0.354 to 0.504 (Table 3). The effective number of codons 
(ENC), which may range from 20 (only one codon is 
used for each amino acid; 1.е., the codon bias is 
maximum) to 61 (all synonymous codons for each amino 
acid are equally used; i.e., no codon bias), were likewise 
intermediate in the primates mitochondrial protein 
condoning genes (Table 3). GC content at the second 
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Table 2 Characteristics for each gene data set in the 54 
Mitogenomic sequences from 44 primate species 


Gene data set N Tiotal Tayn. Tnonsyn. a 





Protein genes 


NDI 954 0.22594 0.56062 0.10797 0.6400 
ND2 1032 0.28488 0.55369 0.19488 0.6829 
СОХІ 1533 0.19501 0.65372 0.04446 0.6678 
COX2 678 0.24309 0.62540 0.11977 0.6908 
ATP8- ATP6 834 0.28099 0.59298 0.17451 0.8017 
COX3 783 0.21378 0.60615 0.08494 0.7019 
ND3 342 0.26373 0.58340 0.15246 0.9531 
ND4L 294 0.25770 0.61033 0.13563 0.6112 
№04 1359 0.25191 0.56225 0.14085 0.6478 
NDS 1788 0.25525 0.56755 0.14766 0.7309 
ND6 486 0.25715 0.50608 0.17487 0.6829 
Cyt b 1134 0.22837 0.57817 0.11319 0.7729 
2. 11217 0.24287 0.58072 0.13166 0.7410 
RNA genes 
12s rRNA 864 0.15497 0.6910 
16s rRNA 1424 0.16676 0.6700 
21 concatenated tRNA 
genes (excluding 1259 0.13533 0.5689 
lysine tRNA) 
All genes 14764 0.22121 0.7229 


N: Numberof sites; mora: nucleotide diversity for all sites; Zyn.: nucleotide 
diversity for synonymous sites; 75,5: nucleotide diversity for nonsynonymous 
sites; o: rate heterogeneity index / Gamma parameter. 


position (СС;) ranged from 38.0% to 40.4%, GC content 
at third codon position (GC,3) ranged from 30.7% to 
50.4%, and overall CG content of the protein codon 
genes (GCc) ranged from 37.4% to 46.5%. Correlation 
analysis indicates that the nucleotide bias affects both the 
synonymous and nonsynonymous sites in protein-coding 
genes. Likewise, there is a highly significant correlation 
between the GC3s and GCc (В? =0.9932, P < 0.001; 
Figure3 A). The correlation between СС; and GCc is 
also highly significant (R?— 0.7772, P < 0.001; Figure ЗВ). 


Neutrality Tests 

As a conservative measure of selection, Z-tests were 
conducted for each protein-coding gene in all species. In 
all cases, dy<dgs was highly significant (P = 0), consistent 
with strong purifying selection of mitochondrial genes in 
these comparisons. The dy/ds ratios differed among 
genes (Figure 4) and these differences may potentially be 
attributable to positive selection or variation in levels of 
functional constraint among genes. The ATP8, ND2, and 
ND6 genes were characterized by а 4/4; ratio well 
above the average, while the COXI gene ratios were the 
lowest among the species. Site-specific positive selection 
inferred from Bayesian posterior probabilities indicated 
that no sites were under positive selection for any gene 
of any of the species. 
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Table 3 Summary of codon usage index of 13 concatenated 


protein genes from 44 Primates species 
Mitogenomic sequences 
Organism ENC  CBI GC, ас;  GCc 


Galago senegalensis 44.147 0.406 0.393 0.390 0.412 


Otolemur crassicaudatus 43.363 0.428 0.399 0.449 0.441 


41.565 0.470 0.397 0.400 0.423 
46.337 0.363 0.394 0.386 0.413 
42.444 0.461 0.390 0.348 0.394 
44.075 0.399 0.383 0.363 0.402 
43.910 0.406 0.399 0.417 0.420 
44.723 0.389 0.389 0.327 0.388 
45.347 0.383 0.394 0.373 0.407 
44.925 0.383 0.393 0.381 0.410 
45.215 0.380 0.394 0.382 0.410 
45.045 0.396 0.391 0.331 0.388 
43.987 0.406 0.390 0.346 0.392 
45.321 0.395 0.392 0.348 0.399 
45.541 0.364 0.388 0.338 0.393 
42.664 0.409 0.383 0.339 0.389 
42.905 0.398 0.385 0.307 0.379 
44.815 0.367 0.385 0.334 0.390 
44.227 0.385 0.385 0.325 0.388 
44.608 0.367 0.387 0.331 0.391 
44.552 0.368 0.387 0.331 0.390 
45.041 0.420 0.401 0.458 0.441 
44.081 0.427 0.397 0.447 0.434 
44.262 0.415 0.403 0.428 0.429 
46.006 0.394 0.382 0.328 0.388 
46.065 0.357 0.389 0.379 0.407 
46.553 0.373 0.387 0.321 0.384 
46.949 0.354 0.387 0.333 0.387 
45.514 0.394 0.384 0.311 0.380 
44.783 0.415 


Nycticebus coucang 

Loris tardigradus 

Perodicticus potto 

Daubentonia madagascariensis 1 
Propithecus verreauxi 

Varecia variegata 

Lemur catta 

Eulemur macaco 

Eulemur mongoz 

Eulemur fulvus mayottensis 
Eulemur fulvus fulvus 

Tarsius syrichta 

Tarsius bancanus 

Saguinus oedipus 

Cebus albifrons 

Aotus lemurinus 

Ateles belzebuth 

Callicebus donacophilus 

Saimiri sciureus 1 

Procolobus badius 

Colobus guereza 

Semnopithecus entellus 
Trachypithecus obscurus 
Presbytis melalophos 

Pygathrix nemaeus 

Pygathrix roxellana 

Nasalis larvatus 

Chlorocebus pygerythrus 1 0.400 0.434 0.431 
Cercopithecus aethiops sabaeus 44.474 0.426 0.399 0.439 0.436 
44.504 0.425 0.396 0.440 0.434 
44.556 0.426 0.400 0.442 0.436 
44.741 0.434 0.400 0.438 0.434 
45.325 0.375 0.387 0.338 0.389 
45.586 0.377 0.380 0.307 0.374 
46.576 0.374 0.387 0.324 0.385 


Papio hamadryas 
Theropithecus gelada 
Macaca sylvanus 
Macaca thibetana 
Macaca mulatta 


Macaca fascicularis 


Hylobates lar 43.281 0.457 0.397 0.498 0.461 
Pongo pygmaeus 43.111 0.436 0.397 0.455 0.442 
Pongo abelii 43.123 0.435 0.398 0.454 0.443 


Gorilla gorilla 1 43.177 0.435 0.398 0.454 0.442 


Pan troglodytes 1 44.396 0.427 0.394 0.456 0.443 


Pan paniscus 41.489 0.504 0.404 0.504 0.465 


Homo sapiens 1 42.622 0.457 0.397 0.474 0.451 


ENC: Effective number of codons (Wright, 1990); CBI: Codon bias index 


(Morton, 1993); GC»: GC content at the second position; GC,3: GC content 


at third codon position; GCc: overall CG content. 
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Figure 2 Sequence diversity of the mitochondrial genome by 

nucleotide position 

Diversity measure was calculated for a sliding window of 500 

bp with a step size of 50 bp. (A) Patterns of total nucleotide 

diversity (mota) for all positions; (B) patterns of 

nonsynonymous nucleotide diversity (Толзул.) for concatenated 

protein genes; (C) patterns of synonymous nucleotide diversity 


(тул) for concatenated protein genes. 


Tests of Recombination 

The measures of linkage disequilibrium are 
summarized in Table 3. The average R (Hill & 
Robertson 1968) across all pairwise comparisons (ZnS) 
for each individual data set and the pooled data provided 
evidence of recombination between relatively distant 
divergent species but not for more recent divergent 
species. The pairwise homoplasy test statistic (Фу) for 
the pooled data also provided similar results to the ZnS- 


Zoological Research 


0.50 





0.45 
са 
5 0.40 
0.35 у=-0.56+2.30х 
Е?=0.9932 
0.30 
0.38 0.40 0.42 0.44 0.46 
GCc 
0.400 
0.395 


0.390 


у=0.30+2.22х 
0.385 Е?=0.7772 





0.38 0.40 0.42 0.44 0.46 
GCc 


Figure 3 Correlation of GC content at the third codon position 
(GC,3), and GC content at the second position (GC;) 
with overall CG content (GCc) 

(A) Correlation of GC,3 and GCc of coding sequences; (B) 


Correlation of СС; and GCc of coding sequences. 


test, but the Фу test statistic provided no evidence of 
recombination for each individual data set. Based on our 
data sets, the assumption of complete linkage of the 
mtDNA genome appears to be justified only between 
more recently divergent species. 


Discussion 


Primate mitochondrial genomes exhibit striking 
patterns of both polymorphism and substitution rate 
heterogeneity among sites and genes, as evidenced by 
sliding-window estimates of sequence diversity, 
nonsynonymous and synonymous divergence, and the 
estimated gamma parameter alpha (о) (Table 2, Figure 2). 
Less primate nucleotide diversity was observed in the 
subunit rRNAs and tRNAs than that in the protein- 
coding genes. Among the 13 coding genes, rates were 
higher in the ND genes and the ATP8- ATP6 gene and 
lower in the COX complex genes, consistent with the 
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Figure 4 Ratios of nonsynonymous/synonymous rate per gene in the mitochondrial genomes of 54 primates. 


Dotted line indicates the average ratio. 


Table 4 Recombination test for primate mitochondrial genomes. 


Internal subclades (node) 








Recombin- 
Gene : 
ation test H, Н, Н; Hy Hs Не 
NDI ZnS 0.0405** 0.0458** 0.0588** 0.0770** 0.2884 0.7266 
Ow 1.2488 1.0927 0.8219 0.6476 0.1496 0.0659 
ND2 ZnS 0.0489** 0.0566** 0.0699** 0.0810* 0.2874 0.7251 
Фу 1.4217 1.2213 0.9230 0.7143 0.1779 0.0579 
COXI ZaS 0.0424** 0.0487** 0.0611** 0.0784** 0.2823 0.7544 
Ow 1.6117 1.4236 1.1128 0.8265 0.1662 0.0345 
cox2 ZnS 0.0795** 0.0874** 0.1147 0.0947* 0.2966 0.7880 
Фу 1.3825 1.1986 0.9053 0.8271 0.1300 0.0000 
ZnS 0.0479%% 0.0558** 0.0696** 0.0845** 0.3012 0.6755 
ATP8-ATP6 

Ow 1.4795 1.2841 0.9917 0.8155 0.1424 0.0772 
COX3 ZnS 0.0407** 0.0553** 0.0575** 0.0730** 0.2273 0.7828 
Фу 1.4295 1.2397 0.9946 0.7820 0.2722 0.0029 
ND3 ZnS 0.0375** 0.0493 ** 0.0610** 0.0790* 0.2893 0.7436 
Фу 1.5716 1.3769 0.9752 0.7419 0.2004 0.0751 
ND4L ZnS 0.0406** 0.0489** 0.0628** 0.0676** 0.2877 0.7621 
Oy 1.3633 1.2171 0.9726 0.7554 0.1937 0.0811 
ND4 ZnS 0.0445** 0.0513** 0.0579** 0.0787* 0.2811 0.7561 
Фу 1.4206 1.2099 0.9199 0.7240 0.1476 0.0101 
NDS ZnS 0.0401** 0.0458** 0.0608** 0.0770** 0.2812 0.7088 
Ow 1.3903 1.2037 0.9311 0.7230* 0.1596 0.0472 
ND6 ZnS 0.0644* 0.0726* 0.0815* 0.0750** 0.2754 0.7926 
Фу 1.264 1.1078 0.9061 0.7237 0.1910 0.0163 
CYTB ZnS 0.0381** 0.0430** 0.0545** 0.0683* 0.2680 0.7586 
Фу 1.3633 1.2238 0.9792 0.7067 0.2383 0.0398% 
AILP. ZnS 0.0442** 0.0504** 0.0640** 0.0770** 0.2835 0.7394 
29; Фу, 1.2291** 1.0311** 0.9620** 0.7483 0.1691 0.0411 


* indicates p«0.05; ** indicates p«0.01 


results from various mammals (Kumazawa et al, 2004). 
Most notably, the nonsynonymous divergence plots were 
similar to those obtained from the pooled data nucleotide 
divergence, while the synonymous changes were only 
slightly different between individual genes of mtDNA 
(Figure 2). Accordingly, rate heterogeneity may result 
from the rate of change at nonsynonymous sites. 

The observed pattern of heterogeneity across the 
mitochondrial genome was similar to the pattern of 
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mitochondrial genome evolution at the intraspecific level 
in humans (Elson et al, 2004) and both the inter- and 
intraspecific levels in Drosophila lard (Ballard, 2000), 
where the synonymous distribution was more 
homogeneous and the nonsynonymous distribution was 
more heterogeneous. Additionally, the results obtained 
here differed from the inter- and intraspecific studies 
based on whole genome in Gadus morhua and its relative 
species (Marshall et al, 2009). Interestingly, all 


Volume 33 Issues E3-4 


Е54 ZHAO, et al. 


mitochondrial genome comparisons mentioned above 
showed significant rate heterogeneity and site specific 
hyper-mutability. Possible explanations for rate 
heterogeneity in the mitochondrial genome include 
selection pressure, variation in the mutation rate of each 
gene (involving base composition bias and properties of 
the L-strand during mtDNA replication) (Marshall et al, 
2009), and the consequence of neutral mutational effects 
(Galtier et al, 2006). The detected evolutionary patterns 
of the mitochondrial genome will influence the choice of 
mtDNA regions used for phylogeographic studies, since 
the more variable regions of the genome harbor excess 
homoplasy relative to the more slowly evolving regions 
(Galtier et al, 2006). Marshall et al (2009) proposed that 
the more variable regions of the mitochondrial genome 
harbor excess homoplasy relative to the more slowly 
evolving regions, so choosing to examine more slowly 
evolving regions for evolutionary analysis among species, 
and even within species, with an increase in the sequence 
length may be wiser. Based on the primate mitochondrial 
genome comparison in this study, we suggest that the 
subunit rRNAs genes, tRNAs genes, and the COX 
complex genes may be promising sites for future 
investigations of primate phylogeography. 

CBIs estimated in this study revealed that there was 
intermediate codon usage bias in primate protein-coding 
genes (Table.3). Several evolutionary processes have 
been postulated as the major factors that determine 
codon usage. For a number of different organisms, it was 
suggested that codon usage was best explained by 
selection for tRNA abundance, gene expression levels, 
and translational optimization (Duret, 2000). In other 
cases, the codon usage was explained by mutation rate, 
mutation preference (Powell & Moriyama, 1997), 
environmental conditions (Goodarzi et al, 2008), 
generation time (Subramanian, 2008), and recombination 
rates (Meunier &  Duret, 2004). Although the 
relationships between codon usage and the underlying 
evolutionary constraints are still not fully understood, 
two major paradigms of codon usage have been found in 
most species (Prat et al, 2009). Natural selection is 
expected to yield a correlation with codon bias in more 
highly expressed genes. Alternatively, codon bias could 
be generated by strictly neutral processes, such as 
mutational biases, where some mutations occur more 
than others across the genome by local variations in the 
base composition (Duret, 2002). Correlation analysis 
indicates that the nucleotide bias affects both the 
synonymous and nonsynonymous sites іп primate 
protein-coding genes. There is a highly significant 
correlation between the GC3s and GCe (В? =0.9932, 
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Figure 3A) and the correlation between GC, and GCc is 
also highly significant (В? =0.7772, Figure ЗВ). Because 
positive selection should not affect silent nucleotide 
positions, our results support the idea that adaptive 
evolution is not the driving force for the codon usage 
bias, indicating that GC mutation pressure may 
determine codon usage in mitochondrial genomes of 
different primate species. 

Using a conservative measure of selection, Z-tests 
indicated that dy«ds was highly significant (P = 0) in all 
cases; this result is consistent with purifying selection of 
mitochondrial genes. A similar result was first found 
during analysis of the mitochondrial genomes of humans, 
chimpanzees, and gorillas (Hasegawa et al, 1998) and 
also reported in additional studies on fruit flies, humans 
and Atlantic cod (Ballard, 2000; Elson et al, 2004; 
Marshall et al, 2009). However, site-specific positive 
selection from a Bayesian perspective indicates that no 
sites are under positive selection for any gene of any of 
the species, indicating that most mutations are nearly 
neutral or slightly deleterious, a pattern consistent with 
the nearly neutral theory of molecular evolution 
proposed by Ohta (1973). Our results therefore indicate 
that the evolutionary pattern of primate mitochondrial 
genomes and each single mitochondrial gene did not 
severely deviate from evolutionary neutrality. 

A presumed absence of recombination is another 
key feature of mtDNA that makes it а useful 
phylogeographic molecule (Marshall et al, 2009). The 
uniparental inheritance of mtDNA indicates that 
recombination should be rare or absent, but a broad 
survey of 267 published animal mtDNA data sets by 
Piganeau et al(2004) provided indirect evidence of 
recombination in more than 14% of cases. Based on the 
pairwise comparisons (ZnS) and the pairwise homoplasy 
test statistic (Dy) for the 13 concatenated protein genes, 
the assumption of complete linkage of the mtDNA 
genome appears to be justified only between more 
recently divergent species. However, no evidence of 
recombination among the individual genes was detected 
with the pairwise homoplasy test statistic (Dw) even 
among the more divergent species, indicating that 
recombination is perhaps a result of substitution rate 
heterogeneity in the primate mtDNA genomes. 
Additionally, as the lineage divergence increases, the 
divergence of each individual gene may produce 
artificial significant ZnS values based on all pairwise 
comparisons from the pooled data and may lead to false 
detection of recombination. Thus, it appears unnecessary 
to consider the effects of recombination on the evolution 
of mtDNA in primates. 
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Description of a new record species of Heterodera from China 
(Tylenchida, Heteroderidae) 


Dong WANG, Lijie CHEN, Yuxi DUAN' 


Nematology Institute of Northern China, Shenyang Agriculture University, Shenyang 110866, China 


Abstract: A new record species Heterodera ripae was found in the rhizosphere soil of nettle (Urtica sp.) on the bank of the Mangniu 
River in Shenyang City, Liaoning Province, China. The new record species can be characterized by vulval cone bifenestrate with 
circular or subcircular semifenestrae, vulval slit length averaging 34.2 um, vulval bridge mostly broad, underbridge weak, bullae 
absent; the second-stage larvae stylet knobs rather wide and slightly projecting anteriorly, tail length averaging 43.7 um and 


hyaline region length averaging 22.3 um. 


Keywords: Heteroderidae; Heterodera; new record species; Heterodera ripae 


INTRODUCTION 


Heterodera ripae Subbotin et al, 2003 (Subbotin et 
al, 1997, 2003), belonging to Nemata, Tylenchida, 
Heteroderidae, was originally found on the roots of 
common nettle (Urtica dioica L.) growing on the bank of 
the Jauza River in the region of Moscow, Russia in 1982 
(Subbotin et al, 1997). Later, the species was also found 
in other nearby European countries, including Estonia, 
Latvia, Armenia, Moldova, Ukraine, Bulgaria, Germany, 
Belgium (Subbotin et al, 1997), Slovakia (Sturhan & 
Lišková, 2004), Greece (Madani et al, 2004), Sweden 
(Andersson & Manduric, 2006) and Spain (López-Robles 
et al, 2011), usually being widely distributed along banks 
of rivers, stream and lakes. During a survey along bank 
of the MangNiu river of Shenyang region in northeast of 
China in 2010, a cyst species was found and identified as 
H. ripae after morphological studies, making it the first 
reported in China as well as outside of the expected 
European environments that had previously been 
identified as habitats. 


MATERIALS AND METHODS 


Populations of the cyst nematode species were 
collected from the rhizosphere soil of nettle (Urtica sp.) 
on the bank of the Mangniu River in Shenyang city, 
Liaoning Province, China. Cysts were isolated from soil 
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samples by the sieving-decanting methods with a coarse 
sieve (openings of 0.18 mm). The second stage juveniles 
and eggs were isolated directly from cysts and eggs were 
hatched. Juveniles were killed by gentle heat, fixed in 
TAF (7 mL 37% formaldehyde, 2 mL triethanolamine 
and 91 mL distilled water) and mounted in anhydrous 
glycerol on permanent slides. Cyst vulval cones were 
dissected from the cysts and mounted in glycerine- 
gelatine. Specimens were examined and measured with 
Motic Images Advanced 3.2, and light microscopes 
equipped with Motic BA400. All measures are presented 
as the mean and the standard error of the mean followed 
by the range in parenthesis. 


RESULTS 


Heterodera ripae Subbotin et al, 2003, new record 
species in China (Figure land 2). 


Measurements 

Cysts (n=10): L=437.6437.4(379.4—488.4) um, W= 
321.7£40.3(262.0-378.6) ит, L/W тайо=1.4=0.1(1.3— 
1.5), fenestral length=47.6+3.6(42.9-53.7) um, mean 
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Figure 1 The illustrations of Heterodera ripae (new record species in China) 
a) Second-stage juvenile head; b) Second-stage juvenile tail; c) The whole of cyst; d) Vulval region. 






Figure 2 Photos of Heterodera ripae by light microscopy 
a) Second-stage juvenile head; b) Second-stage juvenile tail; c) The whole of cyst; d) Vulval region; e) Fenestra and underbridge in vulval region. Black arrow 
indicates underbridge. (a, b, d, e: scale bars=20 um; c: scale bar=200 um). 
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semifenestral width=25.3+2.4 (22.1-29.5) um, vulval slit 
length=34.2+2.3 (31.2-37.6) um, vulval bridge width= 
10.0+1.0(9.0—11.9) um, underbridge length=75.1+3.3 
(70.2—79.7) um, vulva-anal distance=47.145.8(39.4— 
55.2) um. 

Second-stage larvae (n=12): L=400.9430.1(352.8- 
444.9) um, W=20.742.0(18.1-24.2) um, a=19.5+1.9 
(16.7—23.9), b=3.540.3(2.9-3.9), с=9.3+1.6(7.2—12.2), 
с- 3.6-0.2(3.2-4.0), stylet=21.6+0.7(20.2—22.4) um, 
DGO- 4.8+0.4(4.1—5.5) шт, tail=43.7+6.1(33.6—52.2) 
um, hyaline region=22.3+2.6(18.8-26.8) шт, ratio of 
hyaline region to stylet=1.0+0.1(0.9-1.2) . 

Eggs (n=15): L=93.843.1(89.7-100.0) шт, W= 
38.3=1.7(35.7—42.3) um, L/W тайо-2.5--0.1(2.3-2.7). 




































































Description 

Cyst (Figure 1c,d; Figure 2c-e). Lemon-shaped, 
with distinct and rather wide vulval cone, colour varying 
from yellow to pale brown, darkening with age. Surface 
furnished with ridges in irregular zigzag pattern. Neck 
distinct, often forming an angle against body axis. Vulval 
cone  bifenestrate, vulval bridge mostly broad, 
semifenestrae circular or subcircular, vulva set in a 
transverse groove (Figure 1d; Figure 2d,e). Bullae absent, 
underbridge weak, slightly pigmented (Fig. 2e). 
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Second-stage larvae (Figure 1аЪ; Figure 2a,b). 
Body slightly curved ventrally. Labial region rounded, 
about twice as wide as high. Lateral field with four 
evenly spaced lines. Stylet strong, knobs rather wide and 
slightly projecting anteriorly. Median bulb oval, 
pharyngeal glands well developed. Excretory pore 
situated slightly anterior to level of pharyngo-intestinal 
junction, immediately posterior to hemizonid. Tail 
conical, with finely rounded terminus, posterior limit of 
tissues in tail rounded, centrally situated. 

Males. Unknown. 


DISCUSSION 


All morphological data and characters were 
consistent with previous records except a longer J2 body 
on average (400.9 vs 373 um) (Eroshenko et al, 2001; 
Subbotin et al, 1997). 

Heterodera ripae belongs to the Humuli group and 
is closely related to H. humuli and H. vallicola. Though 
it differs from H. humuli in, on average, a shorter J2 tail 
(40-47 уз 49-50 um) and hyaline region (19-23 vs 
26-29 um); and from H. vallicola in, on average, a 
shorter J2 hyaline region (19—23 vs 29 um). 
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